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Abstract 

A derivation of the Hawking effect is given which avoids reference 
to field modes above some cutoff frequency ujc ^ in the free- 

fah frame of the black hole. To avoid reference to arbitrarily high 
frequencies, it is necessary to impose a boundary condition on the 
quantum field in a timelike region near the horizon, rather than on a 
(spacelike) Cauchy surface either outside the horizon or at early times 
before the horizon forms. Due to the nature of the horizon as an 
infinite redshift surface, the correct boundary condition at late times 
outside the horizon cannot be deduced, within the confines of a theory 
that applies only below the cutoff, from initial conditions prior to the 
formation of the hole. A boundary condition is formulated which leads 
to the Hawking effect in a cutoff theory. It is argued that it is possible 
the boundary condition is not satisfied, so that the spectrum of black 
hole radiation may be significantly different from that predicted by 
Hawking, even without the back-reaction near the horizon becoming 
of order unity relative to the curvature. 
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1 Introduction 



The Hawking radiation from a black hole of mass M is most copious at a 
wavelength of order M.Q In this sense it is a long distance effect, whose 
scale is set by the mass of the hole. Thus it is odd that all derivations 
of the Hawking effect refer in some manner to arbitrarily short distances. 
For instance, consider Hawking's original derivation the annihilation 
operator for an outgoing quantum field mode at late times is expressed, via 
the free field equations, in terms of annihilation and creation operators for 
ingoing modes at early times, before the matter has collapsed to form the 
hole. The thermal character of the state at late times is then deduced from 
the boundary condition specifying that the initial state is the vacuum (or 
vacuum plus some excitations of finite total energy.) 

The fishy thing about this derivation is that the frequency of the ingoing 
modes diverges as the time of the corresponding outgoing modes goes to 
infinity. This is because all of the outgoing modes, for all eternity, originate 
as incoming modes that arrive at the hole before the formation of the event 
horizon. An infinite number of oscillations of the incoming modes must thus 
be packed into a finite time interval, so their frequency must diverge. 

Other derivations of the Hawking effect also make reference to arbitrarily 
short distances. A recent derivation by Fredenhagen and Haag is based 
on the form of the singularity in the two-point function as x 

approaches y just outside the horizon. Similarly, arguments based on the 
properties of the correlation functions on the Euclidean continuation of the 
black hole metric assume that the correlation functions have the requistite 
analytic behavior, which involves the form of the short distance singularities. 
Finally, arguments (for conformal fields in two dimensions) based on conser- 
vation of the stress-energy tensor |[ assume the value of the trace anomaly, 
which is the result of regulating a short distance divergence of the theory. 

Since the scale of the process is set by the mass of the hole, it would seem 
that it should be possible to avoid the role of ultra high freqencies much 
higher than in deriving its existence. In a previous paper this issue 
was discussed in detail, and two arguments were offered to support this point 
of view, one involving the response of accelerated particle detectors and one 
involving conservation of the stress-energy tensor. These arguments were 

^We use units with G = c = fi = \. 
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not conclusive but they did make it plausible that the Hawking effect would 
occur even if there were a Planck frequency cutoff in the frame of free-fall 
observers that fall from rest far from the hole. 

It now seems a mistake to focus on a Planck frequency cutoff, since the 
same arguments would support the existence of Hawking radiation as long as 
the high frequency cutoff uJc is much larger than M~^. In the present paper 
it will be shown how Hawking's original analysis can be modified to avoid 
reference to ultra high frequencies. This will require the use of an alternate 
boundary condition, which states roughly that observers falling freely into the 
black hole (starting from rest far away) see no particles at frequencies much 
higher than but less than some cutoff uoc- That this condition implies 
the existence of black hole radiation was implicit in Hawking's original paper 
[H], and was later stressed by Unruh 0. One contribution of the present 
paper is to demonstrate in detail how the derivation can be structured so as 
to entirely avoid invoking the behavior of ultra high frequency modes. This 
analysis involves several sticky technicalities, which we have attempted to 
address as thoroughly as possible. 

This alternate boundary condition is not an initial condition, since it is 
imposed for all times. Moreover, for the reason explained above, it cannot 
be derived from the early time vacuum inital condition. It is in the nature of 
the horizon as an infinite redshift surface that the state of the outgoing field 
modes at low frequencies descends from presently unknown physics at very 
high frequency (in the free-fall frame). Thus the validity of the boundary 
condition cannot be proved within a theory that is only valid below some 
high frequency cutoff. The possibility of justifying the boundary condition on 
energetic grounds will be addressed in section]^. Our conclusion will be that it 
is quite possible the boundary condition is not satisfied, so that the spectrum 
of black hole radiation may be significantly different from that predicted by 
Hawking, even without the back-reaction near the horizon becoming of order 
unity relative to the curvature. Violations of the boundary condition leading 
to a large back-reaction also seem possible, however in such a situation the 
quasi-static, semiclassical framework of our calculations is unjustified. 

The rest of the paper is organized as follows. In section |^ Hawking's 
original derivation is reviewed. In section ^ the role of ultra high frequencies 
in this derivation is discussed, and in section ^ our alternate boundary con- 
dition is formulated and discussed in detail. It is shown in section ^ that this 
boundary condition implies the existence of the usual Hawking radiation. In 
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section ^ the physical basis of the boundary condition is discussed, and the 
imphcations of a violation of the boundary condition are studied. Section 
1^ contains some concluding remarks, and the appendices contain technical 
material needed in the rest of the paper. 



2 Hawking's reasoning 

In this section Hawking's original derivation [||] of black hole radiation from 
a non-rotating, uncharged black hole will be reviewed. We use Wald's formu- 
lation in terms of individual wavepackets, rather than Bogoliubov trans- 
formations between orthonormal bases, because selection of a complete basis 
is distracting and unnecessary for our our purposes. 

Consider an outgoing positive frequency wavepacket P at late times far 
from the black hole, centered on freqency uj and retarded time u. (The 
retarded time coordinate is defined in Appendix A.) Suppose P is normalized 
in the Klein-Gordon norm, so the annihilation operator for this wavepacket 
is given by 

a(P) = (P,$), (1) 

where the bracket notation denotes the Klein-Gordon (KG) inner product. 
(See Appendix B for the definition of the KG inner product, and Appendix 
C for a discussion of this characterization of annihilation and creation op- 
erators.) We are interested in the state of the quantum field "mode" corre- 
sponding to this wavepacket. This is partly^ characterized by the expectation 
value of the number operator, 

{N{P)) = {^\a^{P)a{Pm. (2) 

Using the field equation = 0, this number operator can be expressed in 
terms of operators whose expectation values are fixed by initial conditions or 
other assumptions on the properties of the state |\E'). 

Propagating the wavepacket P backwards in time, it breaks up into a 
"reflected piece" R that scatters off the curvature outside the matter and 
out to past null infinity X~, and a "transmitted" piece T that propagates 



^For simplicity we focus on the expectation value of the number operator. In fact, 
the form of the annihilation operator a(P) discussed below implies also the true thermal 
nature of the state. (See for example ||, 10, 1^.) 
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back through the collapsing matter and then out to X . (See Fig. 1.) The 
original wavepacket P can be exressed as the sum of these two solutions, as 

P = R + T, (3) 

and the annihilation operator for P can thus be decomosed as 

a{P) = a{R) + a{T). (4) 

Since both the wavepackets and the field operator satisfy the wave equation, 
the KG inner products in (P are conserved, and can therefore be evaluated 
on any Cauchy hypersurface. 

Because of time translation invariance in the part of the spacetime ex- 
terior to the matter, the reflected packet R consists of the same frequencies 
with respect to the Schwarzschild time coordinate at X~ as the packet P at 
future null infinity Thus the operator a{R) = {R, $) is an annihilation 
operator for an incoming wavepacket centered on frequency Co. Assuming 
that this mode of the quantum field started out in its ground state, we have 

aiRm = 0, (5) 

so that 

{NiP)) = {^\a\T)aiTm. (6) 

At the transmitted packet T is composed of both positive and negative 
frequency components with respect to the asymptotic Schwarzschild time, 

T = T'+)+T(-', (7) 

and we have the expansion 

a{T) = a{T^+^)-a\T^-^*). (8) 

Thus a{T) is a combination of annihilation and creation operators for incom- 
ing wavepackets at X^. Assuming that both the positive frequency part and 
the complex conjugate of the negative frequency part of the packet T started 
out in their ground states, we have 

a(T'+))|$)=0 a(T(-'*)|*) = 0. (9) 
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Thus, using (^), (Q), (j^) and the commutation relation between annihilation 
and creation operators (49,50) we have 



(iV(P)) 



(T(-',T(-)). 



(10) 



For a wavepacket with a spread of frequencies Acj <^ k, the Klein- Gordon 
norm of the negative frequency packet T'"' can be evaluated in terms of the 
norm of T as described in Appendix D, and one finds, using (57), 



This is just what the emission would be from a body at temperature k,/2tt = 
l/SnM, for a mode of energy O with absorption coefficient {T,T). 

3 Ultra high frequencies 

The difficulty with this analysis is that at past null infinity, the incom- 
ing packet t consists of extremely high frequency components, whose fre- 
quency (with respect to the asymptotic rest frame of the hole) grows as 
~ exp('u/4M)u; as the retarded time u of the outgoing wavepacket goes to 
infinity. This exceeds Planck frequency up for u > AMln{ujp/uj), that is, af- 
ter only several light crossing times for the hole. That the frequency diverges 
in some such manner is immediately evident from inspection of Fig. 1. An 
infinite amount of time at infinity corresponds to the interval between any 
finite u and the horizon at m = oo. The correspondingly infinite number of 
field oscillations must all be packed into the finite range of advanced times 
between some v and Vq, the advanced time of formation of the horizon. 

It is unsatisfactory from a physical point of view to base the prediction 
of black hole evaporation on an assumption that involves the behavior of 
arbitrarily high frequency modes. We are ignorant of what physics might 
look like at those high frequencies or corresponding short distances. In order 
to be confident of the prediction of Hawking radiation, one should formulate a 
derivation that avoids this ignorance while invoking only known physics — or 
at least only more reasonable extrapolations of known physics. 

It is not the unknown physics of high energy interactions that we are 
concerned about here. Although we are dealing with incoming wavepackets 
with arbitrarily high frequency relative to the frame of the collapsing matter 



{T^-\r-^) = {T,T){exp{2nu/^)-l) 



1 



(11) 
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that forms the black hole, there is no interaction between these wavepackets 
and the collapsing matter. The reason is that these incoming wavepacket 
modes are in their ground state, so there is nothing for the collapsing matter 
to interact with. 

What we are concerned about is the need to assume that the physics is 
Lorentz invariant under arbitrarily large boosts. Assuming Lorentz invari- 
ance, one can of course argue that although the frequency of the transmitted 
wavepacket t grows as exp('u/4M) with respect to the asymptotic rest frame 
of the black hole, there is always a local Lorentz frame in which the fre- 
quency appears as low as one wishes. The velocity of this frame relative to 
the black hole approaches the speed of light as -u ^ oo, with a boost factor 
7 = (1 - v2)-i/2 = exp(M/4M). 

We have no observations that confirm Lorentz invariance at the level 
of such arbitrarily high velocity boosts O, O, IT^. Probably the highest 



boost factors at which Lorentz invariance might be checked anytime soon 
arise in cosmic ray proton collisions. We are basically at rest with respect 
to the cosmic microwave background (CMB) radiation. Assuming Lorentz 
invariance, one predicts that for proton energies greater than about 10^° eV 
(relative to the CMB frame), the head-on collision of a proton with a CMB 
photon can produce a pion. This process would leave its mark on the cosmic 
ray proton spectrum. If this mark is eventually observed, it will lend support 
to the assumption of Lorentz invariance that went into the calculation.^ The 
boost factor here relating the CMB frame to the center of mass frame of the 
collision is a "modest" 7 ~ 10^^. 

In the black hole situation, after a retarded time interval Au ~ 4M In 10^^ ^ 
lO^M, the boost factor required to transform an incoming wavepacket to low 
frequency would have increased by more than 10^^. Thus the above deriva- 
tion of a steady flux of Hawking radiation depends on the assumption of 
Lorentz invariance arbitarily far beyond its observationally verified domain 
of validity. 



■^According to Sokolsky it should be possible to confirm this prediction in the 
coming decade. 



7 



4 Cutoff boundary condition 



To avoid the need to make assumptions regarding arbitrarily high frequency 
behavior we will have to give up the attempt to derive the properties of 
the state of the quantum field at late times from the initial condition that 
it is the vacuum state before the hole forms. Instead, we will formulate a 
different "boundary" condition on the state that will still imply the existence 
of Hawking radiation. 

The alternate boundary condition is expressed in terms of the particle 
states defined by free-fall observers near the horizon that have fallen in from 
rest far away from the hole. For frequencies much higher than M~^, these 
particle states are well defined by field modes with positive frequency with 
respect to the proper time of the free-fall observers. Our boundary condition 
will be that outgoing, high freqency field modes are in their ground states. 
How high is "high"? Roughly, to predict Hawking radiation to an accuracy 
?7 <^ 1, it will suffice to assume that the outgoing modes of free-fall frequency 
~ rj^'^M^^ are in their ground state. The statement of the boundary condi- 
tion just given is appropriate for a massless, free field. We defer to subsection 
^4.6| a brief discussion of the modifications required for a treatment of massive 
and/or interacting fields. 

To derive this alternate boundary condition from the condition that the 
initial state is vacuum requires appeal to arbitrarily high frequency modes, 
for the reason discussed earlier. Thus we make no attempt here to derive 
this alternate boundary condition, but rather take it as given. The question 
of physical plausibility of the condition will be taken up in section ^. 

4.1 Precise formulation of the boundary condition 

Actually imposing the alternate boundary condition in terms of the proper 
time of the family of free-fall observers is somewhat complicated. Instead, 
shall employ the affine parameter along radial ingoing null geodesies as the 
relevant "time" variable. This turns out to amount to the same thing near 
the horizon, as will now be explained. 

First note that the usual radial coordinate r is an affine parameter along 
the radial null rays (see Appendix A). To find the rate of change of r with 
respect to the proper time r along the free-fall geodesic, note that the quan- 
tity py = g^^dx^ /dr = (1 — ^)dv/dT — dr/dr is conserved, since the metric 
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is independent of v in Eddington-Finkelstein (EF) coordinates (41). If the 
geodesic starts from rest at oo, one has at infinity dv/dr = 1 and dr/dr = 0, 
so Pt, = 1. It follows then that at the horizon r = 2M, one has dr/dr = —1. 
That is, r is changing at the same rate as the proper time. 

An outgoing solution / to the wave equation near the horizon is nearly 
independent of v in EF coordinates, since the lines of constant r are nearly 
null there. Along the free-fall world line near the horizon, we therefore have 
df/dr = {df/dr)dr/dT = —{df/dr). Thus, for outgoing modes near the 
horizon, the frequency with respect to r on a constant v surface is effectively 
the negative of the frequency with respect to the free-fall observers. 

The particle states of our boundary condition will correspond to wavepack- 
ets / composed of field modes on a constant v null hypersurface S of the form 

fu^imir, 0, 0) = exp{tujr)Yimi9, 0) . (12) 

In an effort to avoid confusion I will call these positive r-frequency modes, 
because they have positive frequency with respect to the proper time of the 
free-fall observers. We can regard the operator a{f) = (/, $) as (proportional 
to) an annihilation operator for a one particle state provided that the Klein- 
Gordon (KG) norm of / is positive. (This is discussed in Appendix C.) 

To evaluate the Klein-Gordon inner product (44) on S, we use the metric 
components in EF coordinates (41) and the surface element (46) to find 
yf—gg^'^d'Li, = —6r'^r'^sin6drd6d(j). Thus the KG inner product takes the 
form 

{f,g) = --jdnj^ drr^rdrg-gdrD- (13) 

This shows that the modes f^^im (0) indeed have positive norm for u; > 0, 
as do localized wavepackets constructed by superposing them.0 

"^It is tempting to try to define a full Hilbert space of one-particle states on a constant 
w-surface using the positive r-frequency modes. However, the fact that the r-integral runs 
only over the interval [0, oo) leads to a problem with this definition. Positive frequency 
modes of the form fujim and fui'im ([l^ ) are not orthogonal for uj ^ uj' , and linear combi- 
nations of positive frequency modes can have negative norm. This is not a problem if one 
restricts attention to wavepackets that have negligible support near r = 0, since for them 
it makes no difference whether the r-integration is over [0,oo) or (— cx),cx)). (One cannot 
take wavepackets of compact support since that would be inconsistent with their being 
composed of purely positive frequencies.) In any case, we will refer to only one wavepacket 
at a time, with no need to consider the full Hilbert space of one particle states. 
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Our alternate boundary condition can thus be implemented as follows. 
We choose to calculate the expectation value of the number operator corre- 
sponding to wavepackets P with the property that on some constant v sur- 
face, V = Vc, their transmitted piece T has only components with r-frequency 
cu^^^ much higher than but less than some cutoff frequency cUc, 

uj, > J"^^ > M-\ (14) 

(If the frequency at infinity uj is much greater than M^^, we also require 
a;'^''^ ^ UJ.) Then, instead of propagating the transmitted piece T of the 
wavepacket P all the way back through the collapsing matter and out to past 
null infinity, we stop when it reaches v = Vc- There we decompose it into its 
positive and negative r-frequency parts and impose the boundary condition 
that the positive r-frequency part (and the complex conjugate of the negative 
frequency part) are in their ground states.^ To carry out this program, it 
must first be established that there exist positive w-frequency wavepackets 
with the property than on some surface v = Vc, their r-frequency components 
satisfy (0). This will be accomplished in subsection ^l3| below. 

4.2 Self-consistency of the boundary condition 

Note that for a wavepacket centered on frequency u and retarded time u, the 
surface v = Vc must necessarily move to the future as u grows with uj fixed, 
in order to avoid the occurence of r-frequency components above the cutoff 
frequency. Thus our boundary condition is not being imposed on a single 
Cauchy surface, so is not an "initial" condition. This raises the question 
whether our boundary condition is consistent with the field dynamics. 

For simplicity, let us think of the boundary condition as being imposed 
on a surface of fixed radius, r = rb.c, just outside the horizon.^ This surface 
is timelike, so the site of the part of the boundary condition imposed at 

^Although the wavepacket P is completely outside the horizon, its positive and negative 
r-frequency parts have support both inside and outside the horizon. (See equations (|30|), 
(53), (54).) 

^Actually, the boundary condition refers to the region inside the horizon as well, since 
the positive and negative frequency parts have support inside the horizon. It is therefore 
more accurate to think of the boundary condition as being imposed on a pair of surfaces 
of constant r, one just outside the horizon and one just inside. Since the one inside is 
spacelike, no question of consistency arises for that part of the boundary condition. 
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advanced time v includes, within its past, sites of parts of the condition 
imposed at earher advanced times. Is the condition imposed at v consistent 
with the earher ones? 

The boundary condition refers to the state of outgoing modes with r- 
frequency to^"^^ in the range cUc > co^'^^ ^ M~^. The modes of frequency 
come from two sources: modes that propagate out from yet closer to the 
horizon with yet higher frequencies, and modes that have scattered off the 
geometry. The state of the former modes can be freely specified, since they 
are above the cutoff until they reach advanced time v and hence no condition 
at all is imposed on them until then. Thus there is enough freedom to con- 
sistently assign the state of the outgoing modes at Uc- But one may still ask 
if the gound state boundary condition is the appropriate one, in view of the 
contributions from the modes that have backscattered. For instance, some 
Hawking radiation can scatter back towards the hole and then scatter again 
out from the hole, apparently leading to some non-zero occuption number in 
an outgoing mode that the boundary condition assigns to its ground state. 
The scattering amplitude for these modes in this region of the spacetime is 
very small however, so such processes should affect the state only very little. 

Now let us consider the modes with frequency less than the cutoff. The 
state of these modes can not really be independently specified, since they can 
be traced back (primarily) to modes yet closer to the horizon with frequency 
Uc, on which a (ground state) boundary condition has already been imposed. 
Thus the state of the modes with frequency u^^'^ < Uc must be calculated, not 
assigned. In fact, it follows from the argument in section ^ that no modes are 
excited while they are propagating close to the horizon; it is not until they 
climb away significantly (on the scale of M) that the presence of Hawking 
radiation becomes apparent in the free-fall frame. 

Thus it appears not inconsistent to impose our ground state boundary 
condition, at least to the order of precision of our calculations. Note that 
we can really only check self-consistency of the calculation: As shown in the 
next two subsections, the unavoidable spread of the wavepackets makes it 
necessary to imose a boundary condition on a wide range of frequencies from 
the beginning. Then all we can do is verify that this boundary condition is 
self- consistent. 
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4.3 Existence of the required wavepackets 

Let p^im denote the solution to the massless scalar wave equation in Schwarzschild 
spacetime that is purely outgoing at future null infinity (and is therefore out- 
going at the horizon as well), and is of the form 

p^i^ = {2nLu)-'/^ exp{-tujt)r-'f^i{r)YUO, 0), (15) 

with 



iLor 




, . I - +A^ie as r* ^ +00 , . 



where r* is the tortoise coordinate defined in eqn. (|4^). These modes are 
normalized according to (p^im,Pt^'/'m') = 5{(jj—uj')5ui5mm.i ■ Using these modes, 
we seek to construct wavepackets that satisfy the condition (p^) restricting 
the r-frequency components on a constant v surface, v = Vc- 
The wavepackets we will employ are of the following form: 



M du B~l exp{iuju) p^im- (17) 

JiU 



Luulm 



Pujulm is a unit norm, positive t-frequency wavepacket centered on frequency 
is a normalization factor, and the factor B^i (inverse of the 
transmission amplitude) is included in the integrand so that we will have 
control over the spread of the part of the packet near the horizon. The 
wavepacket Poiuim is defined by its (purely outgoing) behavior at X"*" and the 
fact that it vanishes on the horizon. Alternatively, propagating it backwards 
in time from X+ as in section ^ one sees that it is generated by data on a 
Cauchy hypersurface formed by a constant v surface v = Vc together with 
the part of X~ that lies to the future of Vc- The wavepacket generated by 
the data at v = Vc alone will be called the "transmitted packet" TQuim, and 
that generated by the data at I" will be called the "refiected packet" Rcjuim- 

Thus we have P^ulm = T^julra + RQUlm- 

For each uj and for u sufficiently long after the collapse that formed the 
black hole, one can always choose Vc sufficiently far in the past so that TquItu 
is concentrated near the horizon. In this case, the asymptotic form f^^i = 
B^iexp{iujr*) can be accurately substituted in the integrand ([T7|) and one 
obtains 



Loulm 



M{2Ti)-^/'^r-^Yim{e, 0) r^^"" duj exp{ioj{u - u)). (18) 
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This transmitted wavepacket is localized in retarded time u, centered 
roughly on -u, with a spread Am ~ Stt/Ao;. More precisely, the spread of Tcjuim 
in u is of course infinite, but the packet is well localized in the following sense.Q 



After carrying out the angular integrals the KG norm (|T^) of the packet (|n 
calculated &t v = Vc reduces to a numerical factor times an integral over 
X of the quantity (sinx/x)^, where x = Auj{u — u)/2. One can show that 
/q (sinx/x)^ dx = (7r/2)[l — (l/7r?/)+0(?/^^)]. Thus, defining rj as the fraction 
of the full norm omitted in a range Am, one has ~ l/iry = A/ttAuj Au, or 

?7~1/AcjAm. (19) 

In Hawking's paper |T|], wavepackets of the form ( ]T7| ) (without the factor 
of B~l) were also employed, however Aco was chosen very small compared 
with the surface gravity k, = 1/4M, so that the wavepackets relevant to 
the black hole radiation would be very peaked in frequency, thus simplifying 
the analysis. From our point of view, the difficulty with this is that such a 
packet cannot be squeezed close enough to the horizon without containing 
r-frequencies above the cutoff Uc- In fact, one must take Auj>k,, and to 
maximize the precision of our derivation one should take Au ~ y^uJj^, as 
will now be shown. 



4.4 Precision of the derivation 

The precision of the derivation we will give is limited by the fact that the 
wavepackets will not be infinitely squeezed up against the horizon. The 
resulting "error" is of order C^ax = (1 ~ SM/rmax), where rmax is the largest 
value of r to occur in the wavepacket .0 Of course, strictly speaking, rmax = 
oo, but a fraction {1 — rj) if the wavepacket is contained within a smaller 
range of r values, given by Au ~ l/rjAu. Thus to minimize the errors we 

^The wavepacket Paaim at 1+ does not have the same width in u as does TQ^im at 
Vc- The wavepacket is somewhat dispersed, since the different frequency components have 
unequal transmission amphtudes. We included the factor B^l in the definition ( p7|) of 
Pouim SO that our packet would be well localized at Vc', it will not bother us that Pouim is 
not as well localized at I"*". 

^Actually, since only a fraction of the wavepacket is located at r ^ ?'max, with the rest 
at smaller values of r, the error is somewhat smaller. To keep the crude analysis that 
follows from getting too complicated, we will simply make the conservative error estimate 
using the largest value of r. 



13 



should minimize the combined error due to the fraction rj of the wavepacket 
beyond r^a.x, and due to Cmax not vanishing. To carry out this minimization 
calculation, we must express Cmax as a function of rj and Au, and minimize 
the error function 

E\v,^uj)^v' + CL.iv,^^)- (20) 

The relation between u and r at constant v is given (cf. (^2]) , (^31) ) by 
du/dr\^ = -2(1-^)-^ = -2C-\ where C = (1-^). It is this factor that 
converts between u-frequency and r-frequency at fixed v, 0;^*"^ = —2C~^uj. 
We assume that on the constant v surface, the wavepacket is squeezed very 
near to the horizon, since that is in any case required in order to deduce the 
existence of Hawking radiation from our boundary condition. Then we have 
(with K = I /AM) 

Cmax/Cmin - exp(ftAM) ~ exp{K/r]Auj) . (21) 

Now assuming the highest r-frequency present in the wavepacket is the cutoff 
frequency, we have = a;^^ = C^iVmax, so that C^in = c^max/^^c- Together 
with (^) this yields 

Cmax ~ exp(fi;/?7Au;) {lj + Au)/uJc , (22) 

where we have returned to the notation u = cjmin- For the purposes of mini- 
mizing the error, we will consider u as fixed, since this is really determined 
by which frequencies we want to learn about. 

Already (^21) shows us that it is not acceptable to choose 6u k a.s 
Hawking did. For instance, suppose that Co ^ 0{k), so the frequencies most 
copious in the Hawking radiation will be included, and suppose that Auj = 
0.01k and rj = 0.01. Then we have Cmax = exp(10,000) k/uJc, which will be 
smaller than unity only if h/uJc is much smaller than we want to assume! 

To minimize the error ( pp] ) , we use ( P^ and set dE /drj = and dE / dAuj = 
0. Up to factors of 0(1), this yields at the minimum: 

7] ~ {Aujf/KUJ^^ and Cmax ~ {Aujf/K^UJ^^ , (23) 

where Aa; satisfies 

cu + Aujc^iAiuf/K^ujl. (24) 
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As long as ^ ^/kU^, the solution is given by 

AUJ ~ y/KuJ^, 7] ~ \/k,/uJc, C'max ~ \[^J^c- (25) 

Note that for such a "minimum error" wavepacket with cj^^^ = cUc, we have 
'^min = ^mLP ~ {u'^ujc/ tc'Y/'^ K , which will Satisfy the condition u^'''^ ^ k as 
long as a; {k/uJc)^/'^ k. 

We conclude that one can work with wavepackets with r-frequencies in 
the required range, with a built-in imprecision of the calculation^ limited to 
an error of order ^JhJuJc. 

4.5 Horizon fluctuations 

Another point that should be checked is how close to the horizon is our 
boundary condition being imposed? If this is within the expected range of 
quantum fluctuations of the horizon itself, then we will not have succeeded in 
formulating a derivation free of short distance uncertainties. To estimate the 
radius rtx. at which the boundary condition is being imposed, note that for 
a mode of frequency coming from a hole of mass M, we have u;*^^'-' ~ ujc 
when (1 - 2M)-i]Vf-i ~ u;^, or rb.c. ~ 2M + Z^, where Ic = uj-\ The scale 
of quantum fluctuations of the horizon 6r can be estimated by using the 
Beckenstein-Hawking entropy S = jA/tl, and setting 55" ~ 1, which is char- 
acteristic of thermal fluctuations about equilibrium.^^ Assuming the horizon 
should be treated as = v4//p independent fluctuating area elements, each 
of area a and radius r, we have 5 A ~ \/N 6a ~ lp6r, so 6 A ~ Ip gives 5r Ip. 
Thus for a Planck scale cutoff, we are perhaps not justified in ignoring the 
quantum fluctuations of the horizon in our derivation. The simple way out 
is to take the cutoff length much longer than the Planck length. This is 
fine until we come to discussing the physical justification for the boundary 

^It may be that the derivation can be improved, reducing the imprecision. The 
wavepacket analysis employed here seems a rather clumsy approach to the problem. The 
problem can also be formulated using the approach of Fredenhagen and HaagQ, which 
focuses on the behavior of the two-point function. That approach may turn out to be 
more suitable for maximizing the precision of the derivation. 

^°This gives the same scale as the one obtained by York using the uncertainty priciple 
and the spectrum of quasinormal modes or using the Euclidean partition function 
approach jlTt . 
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condition, or violations of it. It should be kept in mind that if the modes are 
followed all the way back to where they are squeezed up within one Planck 
length of the horizon, several grains of salt should be added to the whole 
analysis. 



4.6 Massive or interacting fields 

In order to apply the arguments of section ^, it is necessary that the prop- 
agation be governed by the massless wave equation for a sufficiently long 
interval of advanced time v. Thus for a free field of mass m one must impose 
the boundary condition on wavepackets satisfying tu*^^^ ^ m, in addition to 
the condition oj'^'^^ S> already discussed in section ^4.1|. Then one finds 



that particles corresponding to these wavepackets are created near the black 
hole just as are massless ones, and they then propagate away from the hole as 
massive particles. As long as the mass is much less than the cutoff frequency, 
m <^ ujci there is no obstruction to extending our argument to cover the case 
of massive particles. 

It is generally believed that the Hawking effect occurs for interacting fields 
as well as for free fields, although this has never been demonstrated explicitly. 
For the purposes of determining what would be emitted by a real black hole. 



some researchers |]T8| have assumed that the process can be divided into two 
stages, much as for the massive free field just discussed. In the first stage, 
which takes place very near the horizon, the dynamics of the field is governed 
by the asymptotically free regime. In QCD for example, free quarks and 
gluons are assumed to be radiated with a thermal spectrum. In the second 
stage, as the particles climb away from the horizon, the self-interactions of 
the field become important, and the free particle states hadronize into jets. 

A direct demonstration of the validity of this picture has never been 
given, although there are various arguments that support it. Gibbons and 
Perry[^] argued that the periodicity of the Euclidean section of Schwarzschild 
spacetime implies the thermal character of Hawking radiation for interacting 
fields. This argument applies only to the thermal equilibrium state on the 
eternal black hole spacetime. Moreover, it rests heavily on the assumption 
that a state that is regular on the horizon must arise by analytic continu- 
ation from a state that is regular on the (periodically identified) Euclidean 
section. While this condition seems natural in some sense, it has not been 
demonstrated to be necessary. 
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Another argument advanced in favor of thermality is that of Unruh and 



Weiss|]T9|, who demonstrated that the Minkowski vacuum of an interacting 
field theory is a thermal state when viewed by a uniformly accelerating family 
of observers. More precisely, correlation functions in the Rindler wedge are 
given by the thermal density matrix relative to the Hamiltonian that gener- 
ates translations along the boost Kilhng field. This is a purely kinematical 
result. It is, in a sense, a local version of the Euclidean section argument 
that avoids the need for assumptions about regularity of the analytically 
continued correlation functions on the Euclidean section. To turn it into a 
derivation of Hawking radiation for interacting fields, one presumably must 
assume the field is in a state that "looks like" the Minkowski vacuum very 
near the horizon, use the Unruh- Weiss result to describe it from the point of 
view of the static observers as a thermal state, and then propagate this ther- 
mal state out away from the hole. The result will depend on the interactions 
and on what state is incoming from infinity, since this would interact with 
the outgoing Hawking radiation. 

For weakly coupled fields one can study this process using perturbation 
theory. Massless theory in a 2-dimensional black hole spacetime was 
studied by Leahy and Unruh [^], who showed that for an ingoing thermal 
state at the Hawking temperature, the interaction preserves the thermal 
nature of the outgoing state. For an ingoing vacuum state however, the 
outgoing state is not thermal. 

It does not appear to be entirely straightforward to extend the argu- 
ments of our paper to the case of interacting fields, since we use the linearity 
of the field equation to express the annihilation operator corresponding to 
a wavepacket at one time in terms of annihilation and creation operators 
associated with wavepacket at another time. In order to extend our argu- 
ment, one can presumably use the fact that for the first part of the process, 
as the excitations are created, only the propagation of the field "near the 
light cone" is relevant. That is, one can presumably show that only small 
spacetime intervals are involved, and thus use the fact that the correlation 
functions behave like free field ones in this region, due to asymptotic free- 
dom. This picture of the process was outlined by Fredenhagen and Haag in 
the discussion section of 0, but to my knowledge it has never been worked 
out in any detail. 
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5 Hawking radiation in the presence of a cut- 
off 



Having formulated in the previous section a boundary condition on the quan- 
tum state near the horizon that refers only to modes below the cutoff, it is 
now our task to determine the properties of the state far from the hole. 



5.1 Evaluating the occupation numbers 

Suppose now that Poiuim = Ruiuim + TQuim is an outgoing wavepacket of the 
form ([T7|), and propagate TQuim back to a constant v surface v = Vc on 
which its r-frequency components satisfy ujc > o;*-''^ ^ M~^.(More precisely, 
it will be composed of both positive and negative r-frequency modes with 
frequencies in this range.) 

Now we would like to evaluate the expectation value of the number oper- 
ator N{P^uim), subject to our "boundary conditions" on the quantum state. 
These are that 

1. the reflected piece Rouim is in its ground state at X~, and 

2. the positive r-frequency part and the complex conjugate of the negative 
r-frequency part of the transmitted piece T^s/m are in their ground 
states on the surface v = Vc on which the r-frequency components 
satisfy > cu^''^ > 

Subject to these boundary conditions, the evaluation of (A^) goes through as 
in section ^ and we find 

{NiP^uin.)) = -(T^il, (26) 

where T^^lm denotes the negative r-frequency part of the wavepacket T^uim, 
evaluated on the surface v = Vc- 

Now the KG norm in (^6]) cannot have the form of ([TT|) because, as 



explained in sections and the spread of frequencies Au in the packet 
must be taken to be at least of order k, (or even much larger in order to 
maximize the precision). In order to exploit the simple formula (54) that is 
applicable to the negative r-frequency part of a wavepacket of the form (0) 
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with Auj ^ K, we break up the packet Poaim into a large number of pieces, 
defining 

N-l 

Puulra = Pj (27) 

j=0 

l-0+{j+l)Auj/N 

Pj =J\f diu B~l exp(iiu;u) p^im (28) 

Ju}+j/\uj/N 

and the corresponding transmitted packets tj. The {pj} (and the {tj}) are 
an orthogonal (but non-normalized) set of wavepackets, of the type used 
in Hawking' s original derivation when is chosen large enough so that 
Auj/N -C K. (Note that for such large N, B~l does not vary much over 
the range of integration in ( [28| ) and can be pulled out of the integral.) 

Each packet tj has a width of order Am ~ N/rjn, and therefore contains 

r-frequency components in the ratio tu^ix/'^min ~ ^^^^ (see equation (pi])). 
Nevertheless, the full wavepacket T^^im contains only r-frequencies in the 
range (|T1P; the other r-frequency components in the tj's must cancel in the 
sum (|27|), since the sum gives a much more localized wavepacket (which 
suffers much less differential redshift). It is important to stress that although 
we work with the packets tj as a technique to evaluate the r.h.s. of (|26|), we 
do not attribute any direct physical significance or quantum state to them. 
Since extracting the negative frequency part is a linear operation, we have 

/rn(-,r) rp(-,r)\ _ Sr^l,{-,r) J.(-,r)\ /qqN 

V Oulmi -'- wulml — Z^\^ j '''A; /• ) 

To evaluate the KG inner products {tj~'''\t\r''^) we would like to make use 
of the expression (54) for tj"''^' as a linear combination of tj and the "time 
refiected" packet tj. That is, we would like to use the formula 

= c_(e-""^/%- +tj), (30) 

where 

c_ = e-^^^/'=(e-2^'^^/'^ - l)-i (31) 

and 

ujj=to+jALu/N. (32) 

Now this expression for tj"''' was derived in Appendix D assuming that the 
wavepacket tj is squeezed close to the horizon. However, although Tcjuim is 
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squeezed close to the horizon, the individual wavepackets tj may not be, since 
their width Am is much larger than that of TQuim- 

Fortunately this is not a problem, for the following reason. Since the 
KG norm is conserved, we can choose to evaluate (p6D on an earlier surface 
V < Vc, on which not only TQuim but also all the tj are squeezed close to the 
horizon. Moreover, the negative r-frequency part of Tcjuim at w = f c evolves 
to the negative r-frequency part at f < Vc- This is because T^s^m is a function 
of r only through u in this region. Since u = v — 2r — 4Mln(2^ — 1), a shift 
in V is equivalent to a scaling of r near the horizon (where the logarithm is 
dominant) by a linear transformation r ar + b, which leaves the negative 
r-frequency part unchanged. This means we can evaluate the r.h.s. of (^) 
at a surface upon which the tj are sufficiently sqeezed to justify use of the 
formula (pOl). 

The cross-terms in the sum ( P^D vanish, since {tj,tk) = {tj,tk) = {tj,tk) = 
for j 7^ k. The diagonal terms are given by the result (57), so we have finally 

{N{Pa,uim)) = E(^j-.'^.)(exp(27ra;,//€) - (33) 
j 

This is just what the expected occupation number would be for a wavepacket 
mode of the form (^) (equivalently (^)) emitted from a body at temperature 
/t/27r with absorption coefficients {tj,tj) for the component wavepackets pj. 

6 Physics of the boundary condition 

In this section we take up the question of whether there is any way to argue 
that the boundary condition is in fact satisfied. Recall that because of the 
gravitational redshift there is no way, within a cutoff theory, to derive the 
quantum state of the high frequency outgoing modes just outside the horizon. 
The natural expectation would be that they will be in their "free-fall" ground 
state, because from their point of view, there is nothing special about the 
horizon and they are merely propagating along just as they would in flat 
spacetime. The problem with this line of reasoning is that it ignores the very 
question we are trying to address: does the fact that these modes have been 
redshifted down from physics above any cutoff scale leave an imprint on their 
quantum state? 
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6.1 Is this a one-scale problem? 

Together with the presence of the horizon, the absence of any scale other 
than the size of the black hole is really the essence of the Hawking effect. 
One can almost deduce the Hawking result from the fact that the boundary 
condition introduces no length scale other than the Schwarzschild radius 
into the problem. In our form, the boundary condition states that field 
modes near the horizon with r-frequencies satisfying ujc > uJ^"^^ ^ are 
in their ground state. (Since we impose this boundary condition for all 
times, no condition need be imposed on modes with a;*^'^-* > Uc-) This ground 
state is a pure state, however the state of every mode outside the horizon is 
correlated to that of another mode inside the horizon. When only the field 
outside is accessible, there is missing correlation information. An observer far 
from the hole can never determine the state of the modes inside the horizon, 
so the relative phases of the states of all those outgoing modes at infinity 
that emerged from the region of the horizon are completely unknown. The 
state thus cannot be a pure state, but is rather one in which the missing 
information must be maximized in some sense. A maximum entropy state 
is a thermal one, so the state of the outgoing modes should appear thermal 
(modulo absorption coefficients) far from the hole. Since the cutoff Uc plays 
no quantitative role in the problem as formulated, the only scale is M, so 
the temperature must be proportional to 1/M. Calculation shows it to be 
Th = I/SttM. 

In the formulation where the boundary condition is imposed in the asymp- 
totic past, the insensitivity of the black hole radiation to the details of the 
inital state before the hole forms follows from the nature of the horizon as an 
infinite redshift surface: the more time passes, the higher the frequency of 
the relevant ingoing modes. In the limit of infinite time, all that matters is 
the fact that the infinitely high frequency modes are assumed to be initially 
in their ground state. 

But what if one does not assume that physics is invariant under infinite 
blueshifting of scale? If there is new physics at some short distance scale, 
whether it be the Planck scale or something longer, then the gravitational 
redshift may lead to a communication from short to long distance scales 
outside the horizon. That is, the redshift effect leads to a breakdown of the 
usual separation of scales. 

Thus it seems perfectly possible that the quantum state of the outgoing 
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field modes near the horizon might not be the ground state. The precise state 
of these modes could reflect details of physics at much shorter distances. For 
instance, there may be amplitudes for the excited states that could only be 
calculated from a knowledge of the short- distance theory. If this is the case, 
then the spectrum of black hole radiation may be quite different from that 
deduced by Hawking.[^ For example, if one of these modes were to emerge 
at the cutoff in an excited state, then the emission in that mode would be a 
combination of the spontaneous Hawking radiation, the stimulated emission, 
and the original excitation.|^ Thus the flux of energy at infinity would be 
greater than the Hawking flux. 

6.2 Constraints on the stress-energy tensor 

In this subsection we will analyze the implications for the stress-energy tensor 
of a violation of the ground state boundary condition near the horizon. The 
goal is to determine what restrictions energy considerations may place on the 
form of the quantum state of the outgoing modes near the horizon. If the 
components of (T^i/) in the free-fall frame become too large, then neglect of 
the back-reaction is unjustified. I see no reason in principle why this may 
not happen in actuality. It may be that, in fact, the problem of quantum 
fields interacting with gravity in a black hole spacetime defies treatment 
which neglects the back-reaction or which treats it as a small perturbation 
that produces only slow evaporation of the black hole mass. However, if 
this is the case, then the (static) method of analysis used in this paper is 
inapplicable. 

Under what conditions can the back-reaction be treated as a small pertur- 
bation? From the semi-classical Einstein equation Gfj^^ = 87rZ|>(T^i,), we infer 
that the back-reaction will be small provided the stress tensor components 
in the free-fall frame near the horizon are small compared with Ij,^ times the 
typical curvature components there, i.e., 

(T^,) « 1//|M2. (34) 

^^This has nothing to do with the fact that for interacting fields, the spectrum of black 
hole radiation will reflect the dressing and decay of the interacting particle states. Rather, 
we are referring to a difference in the state of the high frequency modes, before the inter- 
actions have had their effect. 

Stimulated emission by black holes is analyzed in Ref. [2^ . 
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In the Unruh or Hartle-Hawking states, one has {T^^) = 0(M~^) in the 
free-fall frame near the horizon, hence in that state the back-reaction is very 
small indeed as long as the hole is much larger than Planck size. In fact, one 
must increase the stress tensor by a factor of order [M/MpY before the back- 
reaction becomes more than a small perturbation. This leaves alot of leeway 
in the form of the state near the horizon, and demonstrates that even within 
the approximation that treats the back-reaction as a small perturbation, 
there is no particular reason why the ground state boundary condition at the 
horizon should hold. 

This leeway in the state at the horizon does not necessarily mean that 
the black hole flux would differ significantly from the Hawking flux however. 
The reason is that the energy carried by outgoing modes near the horizon is 
vastly redshifted by the time they make it out far from the hole. In order to 
make a significant difference in the flux at infinity, an excited outgoing mode 
near the horizon must have a very high energy with respect to the free-fall 
frame. 

To obtain a very crude estimate of the energy density associated with such 
an excited mode, consider a wavepacket that far from the hole is centered 
on a frequency uo with a width Auj ~ uo and a spread in retarded time 
Am ~ uj~^. Suppose this mode is occupied in a one particle state near the 
horizon at some r. As discussed in section |4.3| , its energy relative to the 
free-fall frame will be roughly (1 — ^)~^ uj^ and the proper volume of the 
thin spherical shell containing it will be roughly (1 — — ) uj'^M"^ (since it 
has a thickness Am ~ uo^^ at infinity). Thus the energy density will be 
roughly (1 — ^)~^ c<j^M~^.0 If this mode is followed all the way back to the 
horizon, the energy density diverges, and the neglect of the back-reaction is 
totally unjustified. If on the other hand the mode is foUwed only back to the 
value of r for which the r-frequency is equal to the cutoff cUc, then one has 
(1 — '^)^^ ~ ^cl'^-, and the energy density is 1//^M^. Note that this result 
is independent of cu, even though the extra power emitted cj/Au ~ uP' is not. 

Now if the cutoff represents not just an arbitrary scale beyond which we 
are pleading ignorance, but is rather a physical scale at which the nature of 
propagation might fundamentally change, then it might make sense to halt 

^■^As discussed in section [4.3| , the finite width of the wavepacket leads to a differential 
redshift across the packet, so this simple analysis is too crude to produce reliable numerical 
coefficients. 
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the backward-in-time propagation when the r-frequency reaches uJc- Let us 
entertain this possibihty. 

Suppose then that the energy density near the horizon due to the pres- 
ence of an extra particle in the black hole radiation is given by 1//^M^ as 
suggested by the above computation. More extra particles would just multi- 
ply this by the number of particles, irrespective of their frequency.^ (Note 
however that in order not to over-count the degrees of freedom the inde- 
pendent modes should be spaced in frequency by the spread adopted above, 
AciJ ~ cu.) Similarly, for each particle missing from the Hawking flux, one 
expects a negative contribution to the energy density of the same magnitude. 

Now if the back-reaction is a large effect, then our analysis on the static 
black hole background is actually not correct. We see no reasoning by which 
this scenario can be ruled out, but we can say nothing more about it. If on 
the other hand the back-reaction is small, then at least one of the following 
must be true: 

1. The outgoing modes have only small amplitudes to be not in their 
ground state. 

2. There is near-perfect cancellation between the energy densities due to 
"over-occupied" and "under-occupied" modes. 

3. The cutoff length Ic is much longer than the Planck length. 

It is not even entirely clear that (1) is consistent with a small back- 
reaction, since it only implies a small expectation value for the energy density, 
but still allows fluctuations of order l/PM"^. It would seem to require a full 
quantum theory of gravity to determine whether or not the back-reaction 
could really be neglected in such circumstances. While (2) cannot be ruled 
out, it seems somewhat implausible, since there is no apparent reason for 
such cancellation to occur. Also (3) does not seem very likely, since there is 
currently no evidence of any fundamental length scale other than the Planck 
length. Nevertheless, let us just accept these as the logical possibilities that 
they are. Is there any further difficulty with such a scenario of deviation 
from the Hawking spectrum maintaining small back-reaction? 

^''The estimated energy density breaks down however if the frequency is too low, because 
Au ^ uj~^ will become so broad that the differential redshift across the wavepacket totally 
invalidates the assignment of a particular r-frequency to the packet near the horizon. 
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If the spectrum of radiation is different, but the luminosity is the same 
as the Hawking luminosity, then there must be cancellations of positive and 
negative energy contributions, as mentioned in item (2) above. Although 
this scenario does not seem likely, there seems to be no way to rule it out. 
The possibility that the net luminosity differs from the Hawking luminosity 
appears to be somewhat constrained however by general properties of the 
stress energy tensor if the back-reaction is to remain small. 

As first shown in the 70's [§, |^, given some relatively "theory-independent" 
constraints on the behavior of the stress-energy tensor one can derive a for- 
mula for the net radiation flux far from a (quasi)static black hole. These 
constraints are: 

• {T,u)'^ = 0; 

• (T^j/) is nonsingular on and outside the horizon (in regular coordinates); 

• (T^jy) is static and spherically symmetric (in four dimensions); 

• no radiation is incoming from infinity at late times. 

If all these properties hold then it can be shown that the luminosity L of 
the black hole is given in two spacetime dimensions by 

1 roo 

L = -M drr-^T:) {D = 2) (35) 
and in four dimensions by 

/•OO POO 

L = 2-kM dr (T^) + 47r / dr {r - 3M) {T^g) {D = A). (36) 

Let us consider first the two-dimensional case. Then the luminosity is 
determined entirely by the trace of the stress-energy tensor. If we consider a 
conformally invariant massless scalar field, the trace is determined in a state- 
independent manner by the trace anomaly to be (T") = i?/247r where R is the 
Ricci scalar. Putting this in (^) yields the Hawking flux Lh = l/7687rM^. 

Any deviation from the Hawking flux for a conformally invariant field in 
two dimensions thus implies that at least one of the properties of the stress 
tensor assumed above must fail to hold. It seems that the most questionable 
assumption is that of the value of the trace. But what would be the physical 
basis for a deviation from the usual trace anomaly formula? 
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It was argued in Q that the presence of a high frequency cutoff Uc is only 
likely to affect the value of the trace by terms of order 0{R/uj1). This argu- 
ment was based on the assumption that the origin of the quantum violation 
of conformal invariance can be located entirely in the regulated functional 
measure in the manner of Fujikawa |2^. If correct this implies that the cor- 
rections to the trace (and to the Hawking flux) are very small indeed for holes 
much larger than the cutoff length. However, if there is fundamentally new 
physics at the cutoff scale, then the violation of conformal invariance will 
not be due simply to the non-invariance of the regulated functional measure. 
This opens up the possibility of a more significant deviation from the usual 
trace anomaly. Nevertheless, the fact that the usual trace anomaly is state- 
independent (assuming the state has the usual short distance form down to 
some cutoff much smaller than the radius of curvature of the spacetime) sug- 
gests strongly that no significant deviation from the usual trace would occur. 
Thus, at least in this two-dimensional model, it is hard to see how the flux 
could differ from the Hawking flux and still have a small back-reaction. 

In the four dimensional case ( p6D the situation is perhaps different. For a 
conformally invariant field the trace is still determined by the trace anomaly, 
and is given by 

(T„") = /5CV48 = /3MVr^ (37) 

where eOvr^/? = 1, |, SS/eOTr^ for fields of spin-0, 1/2 and 1 respectively, and 
is the square of the Weyl tensor. Now however, the trace does not suffice to 
determine the luminosity. One free function of r remains undetermined. The 
reason is that, unlike in two dimensions where all metrics are conformally fiat, 
the Schwarzschild spacetime is not conformally fiat, so even a conformally 
coupled field scatters in a non-trivial way. Both the spin of the field and 
the detailed radial dependence of the metric affect the radial dependence of 
{T^y) and the net flux at infinity. Numerical computations p5| show, for 
example, that for a massless, minimally coupled scalar field in the Unruh 
vacuum in Schwarzschild spacetime, the contribution of the second integral 
to (^) is relatively small, and the Hawking luminosity is of order Lh ~ 
(48007rM2)-i ~ lO-^M-2. 

A deviation from the Hawking luminosity could be produced, as in the 
two dimensional case, by a deviation from the usual trace anomaly, however 
the same arguments as given in that case make this seem unlikely. But in 
four dimensions there is another possibility: Any change in the tangential 
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stress (Tq) will entail a change in the luminosity L, without violating the 
above assumptions on the behavior of {T^v)- Can this be exploited to allow 
for a deviation from the Hawking luminosity? While it is not clear why there 
should be any fundamental difference between the two and four dimensional 
cases with regard to the possibility of deviating from the Hawking radiation, 
let us just take the result (|36D and see what can be done with it. 

Note first that the r-dependence of (T^ ) has alot to do with the scat- 
tering behavior of fields propagating in the Schwarzschild geometry. Thus 
at most, we should think of the possibility of freely modifying (Tg ) at one 
point, letting the behavior everywhere else be determined by the scattering 
off the background geometry. A change in the luminosity of order 5L could 
be produced, consistent with (P^D, in two qualitatively different ways: (i) a 
change 5Tq ~ 0{SL/M~'^) over a range 6r ~ M, or (ii) a very large change 
6Tg over a very small range of r near the horizon. The second way seems 
inconsistent with the scattering behavior of the field, since the effective po- 
tential that governs the scattering is well behaved near the horizon. The 
first way requires only a relative change 6Tg/Tg of order unity to change the 
luminosity by order unity. Thus there seems to be no obstacle to the physics 
at the cutoff scale leading to a deviation from the Hawking luminosity, even 
if the back-reaction is to remain small. 

7 Conclusion 

What has been accomplished in this paper? We have succeeded in formulat- 
ing a derivation of the Hawking effect (for massless free fields) that avoids 
reference to field modes above some cutoff frequency in the frame of the 
free-fall observers that are asymptotically at rest. To stay below the cutoff 
it is necessary to impose a boundary condition on the field near the hori- 
zon for all times. The boundary condition states roughly that the outgoing 
high frequency field modes are in their "ground state" as viewed by free-fall 
observers. This boundary condition is not derivable from the initial state 
within the cutoff theory. 

The precision of our derivation is controlled by the ratio of the cutoff 
length to the Schwarzschild radius of the black hole, and is limited by \JljM. 
For a black hole large compared with the cutoff length, the largest source of 
imprecision is the unavoidable spread of the wavepackets employed, and the 
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associated large differential in the redsliift suffered across the packet when it 
is near the horizon. 

The boundary condition we impose may or may not be physically the 
correct one. If it fails to hold, then there will be a deviation from the Hawking 
spectrum. It seems this could occur either with or without a large back- 
reaction. If the back-reaction is to remain small, then either the deviation 
must be small, or there must be cancellation between positive and negative 
energy contributions, or there must be a physical cutoff much longer than the 
Planck length. In four dimensions, the generally expected behavior of the 
stress-tensor cannot be used to definitively rule out any of these scenarios. 

Even a very small deviation from the thermal nature of the Hawking 
radiation would seem to entail a breakdown in the generalized second law of 
thermodynamics ^ . Thus one has reason to suspect that the physics 



at the cutoff scale somehow conspires to produce precisely the "thermal" 
state. However, that is not to say that the ordinary effects of quantum field 
propagation in the black hole background should not leave their mark on the 
radiation. The scattering of wavepackets by the geometry is one well-known 
aspect of this mark, but it is conceivable that the redshifting of the physics 
at the cutoff is another one. If that is the case, then the thermodynamic 
behavior of physics in a black hole spacetime may turn out to be much more 
subtle than was previously thought. 

Given a candidate theory with a short distance cutoff, it will certainly 
be interesting to study its behavior in a black hole spacetime, in which the 
redshift effect acts as a microscope to reveal consequences of short-distance 
physics at larger scales. 
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Appendix A: Black hole line element 



The static, spherically symmetric black hole hne element in Schwarzschild, 
tortoise, double-null, and ingoing Eddington-Finkelstein coordinates takes 
the following forms: 

ds^ = (l_2M)di2_(i_2M)-idr2-rW (38) 

= (l-^)(di'-dr*')-rW (39) 

= (l-^)dudv-r^dn^ (40) 

= (l-^)dv^-2dvdr-r^dn^ (41) 

with 

r* = r + 2Mln{^-l) (42) 

u = t-r*, v = t + r*. (43) 

The coordinates u and v are called the retarded and advanced time coordinates 
respectively. 

If x'^(A) is an affinely parametrized geodesic, then it is a stationary point 
of the integral / gfij.x'^x^dX, where the dot • = d/d\. To see that r is an 
affine parameter along ingoing radial null geodesies, it is convenient to use 
the Eddington-Finkelstein coordinates, so that v = 9 = (j) = {). Upon varying 
f (A), one immediately finds f = 0, so r = aA -|- 6 for some constants a and h. 



Appendix B: Klein-Gordon inner product 

The Klein-Gordon inner product (/, g) between two initial data sets / and g 
on a Cauchy surface E is defined by 

{f,g) = lfdE^ (44) 

f = lV^9''{rd.g-gdJ*). (45) 
The surface element dT,^ is given by 

dE^ = ^e/iijfc da'-da^da^, (46) 
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where a' {i = 1,2,3) are coordinates on the surface S. For solutions of the 
KG equation of compact support, (|^) is independent of the Cauchy surface 
on which the integral is evaluated, since the current vector density j'^ is 
divergence free, d^j'^ = 0. We shall have occasion to evaluate the KG inner 
product on a surface that is null, which can be thought of as a limiting case 
of Cauchy surfaces. 

Appendix C: Quantum field theory 

The field operator $ for a real, free scalar field is a Hermitian operator that 
satisfies the wave equation = 0. We define an annihilation operator 
corresponding to an initial data set / on a surface S by 

a(/) = (/,<f)E. (47) 

If the data / is extended to a solution of the wave equation then we can 
evaluate the KG product in (|47|) on whichever surface we wish. The hermitian 
adjoint of a{f) is called the creation operator for / and it is given by 

at(/) = -(/*, (48) 

The commutation relations between these operators follow from the canon- 
ical commutation relations satisfied by the field operator. The latter are 
equivalent to 

Hf),a\g)] = {f,g), (49) 

provided this holds for all choices of / and g. Now it is clear that only if / 
has positive, unit KG norm are the appelations "annihilation" and "creation" 
appropriate for these operators. From ( ^9] ) and the definition of the KG inner 
product it follows identically that we also have the commutation relations 

Hf),aig)] = -{f,g*), [a\f),a\g)] = ~{f\g). (50) 

A Hilbert space of "one-particle states" can be defined by choosing a 
decomposition of the space S of complex initial data sets (or solutions to 
the wave equation) into a direct sum of the form S = Sp ® Sp*, where all 
the data sets in Sp have positive KG norm and the space Sp is orthogonal 
to its conjugate 5*^*. Then all of the annihilation operators for elements of 
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Sp commute with each other, as do the creation operators. A "vacuum" 
state corresponding to Sp is defined by the condition a(/)|\E') = for 
all / in Sp, and a Fock space of multiparticle states is built up by repeated 
application of the creation operators to 

(Instead of thinking of the Hilbert space as the Fock space corresponding 
to some decomposition SpQ)Sp* as above, it is perhaps conceptually preferable 
to take the point of view of the algebraic approach to quantum field theory 



2^, according to which a "state" is simply a positive linear functional p 
on the T^r-algebra of field operators. Thus for example, to express the idea 
that a given field mode / is in its ground state, one says that the state 
p satisfies p{Oa{f)) = for all operators O. This language is preferable 
if, as is often the case for quantum fields in curved space, one wishes to 
simultaneously consider a state as an element of two completely differently 
constructed (for example "in" and "out") Fock spaces. In the algebraic 
approach, no mysterious "identification" of the two Fock spaces is required. 
Another advantage is that whereas the statement that the field operator 
is "hermitian" is meaningless until the Hilbert space on which it acts has 
been specified, the statement that $ goes into itself under the abstract -k 
operation is always well defined. The algebraic approach is clearly preferable 
in contexts (e.g. |2^, ^]) in which one wishes to obtain results valid for a 
class of quantum states that is as wide as possible.) 



Appendix D: Negative frequency part 
of the transmitted wavepacket 

Consider the transmitted part Iqu of a wavepacket p^z^u propagating in the 
Schwarzschild black hole spacetime, narrowly peaked in w-frequency about lD 
(at large r) and about some late retarded time u. For the original Hawking 
argument one needs to determine the KG norm of the negative frequency 
part of tcju at I~ in terms of the norm of t^u itself. For our argument in this 
paper, it is the negative r-frequency part on a constant v surface that is of 
interest. It was Hawking's original argument that these two are related, using 
the geometrical optics approximation to propagate the very high frequency 
modes in question back out to . 

Consider a collection of null surfaces, wavefronts for such a mode. In Fig. 
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2 one surface is shown that is outgoing at retarded time u and ingoing at 
advanced time v{u). The key fact is that for late retarded times, the value 
of the affine parameter r(M,fc) where this wavefront intersects the surface 
f = f c is linearly related to the advanced time f Therefore the negative 
r-frequency part of t^^ aXv = Vc propagates back to the negative w-frequency 
part at I~ . Thus the corresponding KG norms are identical, so in both cases 
we can carry out the calculation ai v = Vc- 

Now there is an observation p, Q that makes the extraction of the nega- 
tive frequency part simple: let U be defined by km = — ln(— /tf/), and consider 
the functions q and q, defined by 

rrn / e""" forf/<0 

= I for t/ > (^1) 

and 

qiU) = q{-U). (52) 

That is, q is just the function q reflected over the hne t/ = (m = oo). Then 
one can easily show that the functions 

= c+(g + e""'"/''g) and (53) 
g'-^ = c_(e^"'^/'^g + g) (54) 

are pure positive and negative [/-frequency packets respectively. One can 
solve for the normalization factors c+ and c_ by setting q = g*"""' -|- q^~K This 
yields 

c_ = -e-'^-/'^c+, (55) 
c+ = (1 - e-'"'^/")-\ (56) 

Finally, the KG norm of g'"' is calculated from (54) and (55,56) using (g, q) = 
— {q, q) and (g, g) = 0, yielding 

(g(-),g(-)) = -(g,g)(e2W«_l)-l. (57) 



^^Hawking argued that this is because as one goes from (u, Vc) back along the wavefromt 
and out to I~ , the "vector" that connects the wavefront to the horizon (and earher to 
the null ray that becomes the generator of the horizon) is parallel transported into itself. 
This is not actually correct, since the connecting vector satisfies not the parallel transport 
equation but the geodesic deviation equation. Nevertheless, one still obtains a finite linear 
scaling of the connecting vector, which is all that is required for the argument p9[. 
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The preceeding calculation is directly applicable to the wavepacket t^^, 
squeezed near the horizon on the surface v = Vc- The relation between r and 
u along V — Vcis given by (42,43), u — Vc — 2r* — Vc — 2r — 4Mln(2^ — 1). 
For our wavepacket near the horizon, the spread in r is very small compared 
with 2M, so the wavepacket only has support where one has km ~ — ^'n.ij^ — 
1) + const.. Thus U and r are linearly related (via —kU = — 1), so 
the negative r-frequency part g*" '^' is equal to the negative ^/-frequency part 
g'"-* (54) with uj — uj, provided the packet is sufficiently peaked in frequency 
about oj (Ao; <^ k) so that the expressions (53,54) for the positive and 
negative frequency parts still hold. This is the result used in the text. 
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FIGURE CAPTIONS 



Figure 1. Conformal diagram depicting wavefronts of the wavepacket P — 
R + T propagating in the spacetime of a spherically symmetric collapsing 
body. 

Figure 2. Conformal diagram depicting the propagation of a wavefront. The 
point (fc, r{u, Vc)) is connected by a radial null geodesic to a point on I~ at 
advanced time v{u). The afhne parameter r along the line v = Vc is linearly 
related to v{u) for late retarded times u. 
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